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ABSTRACT 



A computer-controlled system for automated, remote 
measurement of phase and amplitude shaded doublet, triplet, 
and quadruplet acoustic arrays beam patterns was developed 
and tested. Experimental data were collected with an 
anechoic chamber to minimize backscatter and other 
interference effects. 

Preliminary impedance, reactance and resistance plots 
were made on several speakers to ensure that a matched set 
with similar characteristics was utilized for the arrays. 
Once the arrays were established, a computer program was 
developed for an HP-86 computer which allowed it to rotate 
the arrays, sample and collect data, and plot the results. 

Various combinations of amplitude and phase shading were 
applied to the dipole, tripole, linear quadrupole and 
quadratic quadrupole array. The results were then compared 
to computer generated beam patterns which agreed to i 2dB. 
This uncertainty was mainly generated by the slight 
deviations in the positioning of elements. 
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I. INTRODUCTION 



A. PURPOSE OF STUDY 

The purpose of this investigation was to develop and 
verify a computer-controlled experimental system to measure 
and plot the pressure fields of various acoustic arrays in 
an anechoic environment. This was to be done using various 
combinations of phase shading and amplitude shading of 
doublets, triplets, and quadruplets. Experimentation was 
conducted under computer control to maximize the amount of 
data available and to ensure the highest degree of accuracy. 
Additionally, a comparison was conducted with the 
theoretical predictions obtained by Park [Ref. 1]. 

B. BACKGROUND 

The use of dipoles and quadrupoles for beam forming and 

direction finding, both in air and in water, has been a 

subject of intense investigation since the 1940's. The 

basic radio direction finding apparatus of the 1940's was 

nothing more than a set of orthogonal dipoles which 

consisted of two crossed loops. The relative signal 

amplitude from the two dipoles was used to determine the 

direction of arrival of a signal. The angle is equal to 

Arctan (y/x) where x was the amplitude of one dipole and y, 

c 

the amplitude of the other dipole. This presented a 180 
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bearing ambiguity which was resolved by comparing the phase 
of the X and y dipole signals with the phase of a co-located 
omni-directional antenna. [Ref. 2] 

In the 1950's, extensive studies were conducted using 
planar arrays with multiple acoustic and electromagnetic 
dipole elements. These studies led to the formulation of 
planar arrays employed in radar systems of the 1960's. 
Further research resulted in the ability to steer beam 
patterns by use of multiple elements, phase shading, and 
amplitude shading. 

In 1976, work was done by Hughes and Thompson to 
investigate acoustic communications, particularly in 
underwater acoustics for cases where the position of the 
source or the receiver is fixed. They investigated the 
tilted directional response patterns below the horizontal 
plane due to amplitude weighting and a single 90° phase 
shift. Here, a method of achieving a tilted pattern from a 
symmetrically shaped array without using numerous phase 
shifts or time delays was discussed. [Ref. 3] Moses and 
Smith also examined beamforming this time with a cylindrical 
array using real (rather than complex) shading coefficients 
[Ref. 4]. 

In 1985, Flanagan examined the use of computer-steered 

microphone arrays which were two-dimensional planar arrays 

of 9 X 7 elements. Here an array was designed and computer 

steered to discriminate against sound arrivals from all but 

the desired directions. [Ref. 5] 
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One example of the Navy's current use of dipole 
technology is in sonobuoys. The SSQ-53 DIFAR Sonobuoy 
possesses both an omnidirectional and bi-directional 
hydrophone capability. The SSQ-77A VLAD Sonobuoy possesses 
an upper and lower array strung in a single line array. 
Beamforming is accomplished by the use of both phase shading 
and amplitude shading. 

These advancements in the study of arrays ranging from a 
single dipole to large scale planar arrays, coupled with 
computer control, form the basis for this further study of 
shaded arrays. There is currently a related research effort 
ongoing in the study of computer controlled arrays and beam 
patterns being conducted by the Environmental Physics group 
at the Naval Postgraduate School. 

C. REPORT OVERVIEW 

In Chapter II, a basic introduction to theoretical 
pressure and directivity calculations of a dipole is 
developed. This is then used as a base to expand to the 
physical case and the electrical linear systems case. 
Finally, these results are applied to different forms of 
quadruplets and triplets to show theoretically, the 
expected results. 

Chapter III contains a summary of the beam pattern 
measuring system, the computer program development, 
experimental procedure, and representative data and plots. 
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Additionally, it provides a description of the system which 



was assembled 


for measuring 


the pressure 


fields of 


the 


various arrays. 










Chapter 


IV 


contains a summary of data 


analysis 


and 


ver if ication 


• 


Additionally, 


representative plots 


are 


provided for 


amplification of 


arguments . 






Chapter 


V 


contains 


concluding 


remarks 


and 



recommendations for further work. 

Appendix A provides a brief description of each piece of 
equipment and summary of its capabilities as applied to this 
experimental investigation. 

Appendix B lists the BASIC computer language program 
used for data acquisition and formulation in this study. 

Appendix C contains specific plots, both experimental 
and theoretical, obtained as a result of this investigation 
and that of Park [Ref. 1]. These plots were obtained for 
beam patterns of various array configurations. 
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II. THEORY 



The dipole has been extensively explored and is now 
utilized as a basic unit in the construction of a great many 
arrays and antennae. In this section we will elucidate the 
theory using two approaches 1) an acoustic approach in 2 
dimensions and 2) a linear systems approach in three 
dimensions. These two approaches are related to each other 
and yield the same result after minor changes are made in 
the geometry of the problem. Both are presented because of 
direct Navy applications dealing deal with SONAR arrays. 

Following the presentation of each of these, the 
equations for the triplet and quadruplet arrays utilized in 
this research effort will be shown. 

A. DIPOLE RADIATION-ACOUSTIC APPROACH 

Because it is usually desirable and useful to be able to 
determine the direction from which an acoustic signal 
originated, arrays were developed to provide this 
information; a simple example of this is the acoustic 
dipole. A simple dipole consists of two simple sources of 
equal strength but opposite phase, separated by a distance d 
Fig. (2.1). 
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Figure 2.1 Geometry used in Deriving the 

Radiation Characteristics of an 
Acoustic Dipole. 



The pressure, p(r,©,t), measured at the field point (r,©) 

due to source 1 and 2 respectively are 
Pj = [A;(r 4 - Afj)] exp[j((Ot - k(r +- Afj)] 

P 2 = [-A'(r - Ar,)] exp[j(o)t - k(r + Ar,)] (2) 

where A is the pressure amplitude of an individual source 
measured at one meter from the source, r is the distance 
from the field point to the midpoint between the sources, 
0)“2 77 f is the angular frequency and k is the wave number. 

The minus sign associated with P2 accounts for the 180 
phase difference. The total acoustic pressure is the sum of 
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P| and 



P(r.0,t)r(A'r){([exp(-jkAr,)] [lt(Arj 'r)J) - {[exp(-jkAr 2 )]'[l* (Ar,'r)]j]exp(j(o)t- kr)) ^3) 
Considering only field measurements where the separation d is 

small compared to the distance r, we can use the far field 
approximations, ^'' 1,2 ~ (<J/2) sin 0 and Atj ,/r < < 1 

This yields the acoustic dipole radiation pressure field, 

P(r,0.t) = -j (2A/r) sin {(kd 2)sin 6)} exp j(o>t - kr) 




(a) 




Figure 2.2 Cartesian and Polar Representation 
of Pressure Amplitude Distribution 
of Dipole Radiation. [Ref. 6] 

17 



As noted in Eq. (4) and seen in Fig. (2.2), plotting the far 
field for d<<Ar pressure nulls exist when (kd/2) sin© =77>^ 
where n is an integer. (Ref. 6] 

The twin lobes and nulls of a dipole result in an 18^* 
ambiguity when trying to localize a source. If the dipole 
is combined with an omnidirectional receiver, the ambiguity 
can be resolved. In Fig. (2.3), various combinations of 
dipole and omnidirectional receivers are combined to provide 
a cardioid response. 




Figure 2.3 Directional Diagrams of Various Combinations 
of Bidirectional and Nondirectional 
Microphones. [Ref. 7] 

Here is the voltage output of the omnidirectional 

receiver and R 2 the voltage output of the directional dipole 
receiver. It is not important for the two sources to be of 
the same sensitivity, but the ratio R^/Rxmust be unity for 
true cardioid response. The phase shift between these two 
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receivers, which yields directionality, is achieved through 
a difference in the phase characteristics of the omni which 
is differentiated (high pass filtered) to compensate for the 
linear frequency dependence of the dipole sensitivity. [Ref. 
7] . 

This design is presently used in SSQ-77A VLAD Sonobuoys 
where a cardioid is formed pointing in the downward 
direction. This maintains the null in the upward direction 
in order to alleviate high frequency interference from wind 
and sea surface noise. 

B. LINEAR SYSTEMS 

Consider the geometry for a dipole as shown in Fig. 
(2.4) ; 




Figure 2.4 Dipole Geometry [Ref. 8] 
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In this representation as developed in Ref. (7) , Cj and C_]^ 
are complex weighted point sources which are free to rotate 
in the x-y plane and d is the distance between elements. 

The complex weighted response factor ci point 

source is 

c = 3i exp(j(p.) 

where ai is the amplitude and ^i phase angle. For a 

dipole, the complex weighted response factors are Cj = aj 

and “ -3j « 

For a wave originating from a distant source and 
propagating to the receiver, it is observed that there is an 
essentially constant pressure magnitude over the dipole 
aperture. Because the phase varies with position along the 
wave front in a traveling wave, the response of the dipole 
is given as a function of frequency and direction cosine u 
where u = sin©cosNk* The response is 

D(f,u) = Cj exp(jndu'>.) + c, exp(-jTtdii X) (6) 

This equation is simplified for the dipole to; 

D(f,u) = 2ja, sin (Ttdu/X) (7) 

Since u is a direction cosine, it is only valid over the 
interval -Ku<l. As a result, 0<0<7t and . u = o 

implies a broadside or on-axis geometry and u = 1 implies 
an end-fire geometry. A plot of D(f,u) as a function of u 
can be used to approximate the response of a dipole. [Ref. 
8 ] 
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An example of this two element response of two point 
sources with a separation d equal to one half the wavelength 
of the traveling wave results in a real response for values 
of u from (-X;2d) to ( + X/2d) • If the sources are assigned 

dipole weightings, substitution of the functions © and for 
direction cosine u result in the unnormalized, far field 
beam pattern 

D(f,0.vi/) = j2a, sin[(Kd/X) sin 0 cos \\i] (8) 

This linear systems theory can be easily applied to the 
axial quadrupole and Tesseral quadrupole. The axial 
quadrupole is an array of three complex weighted point 

sources where c, = aj = c.j = a., = 1 and = a^ = -2 as in 

Fig. (2.5). 

Here can also be thought of as 2 sources with each having 
an amplitude weight equal to -1. 




Figure 2.5 Axial Quadrupole Geometry [Ref. 8] 
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with the 1,-2, 1 amplitude weighting, the far-field, 

unnormalized directivity function is defined as: 

D(r,f^) = [exp(-j2n:f^d 2)] - 2 + [exp(j2;:r^d 2)] ( 9 ) 

where = u X = sin6cos<[/ X . This reduces to: 

D(f.0.v|/) = 2 [cos((rtd;X) sinOcosy) - 1] (10) 

In the case of a Tesseral quadrupole Fig. (2.6), a 
planar array of four amplitude-weighted point sources, the 
unnormalized directivity function is: 

= 2Q cos(;rd(r^ + fy)] - 2Q cosl7rd(f^ - fy)] (11) 

where = u X = sinGcosy'X ^ Ty = v/X = sinGsinvt/ X and the 

direction cosine v = sin © sin'^. This yields 

D(f,G,\|/) = -4Q sin[(Tid'X) sinGcosv;/] sin[(nd X) sinGsiny] (12) 




Figure 2.6 Tesseral Quadrupole Geometry 
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which is the far field, unnormalized directivity function 
for this quadrupole as a function of© and^. [Ref. 8] The 
above development is particularly important because this 
type array was investigated in the course of this work. 

C. GENERAL QUADRUPOLE 

Park, in his thesis (Ref. 1), developed an equation for 
calculation of a four element array Fig. (2.7). The 

equation was developed by adding complex pressures from four 
sources to yield 

p = ei(wt - kR - kDjSinSsintP) ^ pj('*^ - kR- + kD^sinGsinO - (p.,) ( 13 ^ 

R R - ^ I 

4 - ^ kD,sin0cos<I> - (p,) 

R 

+ * kD^sin0cos<I> - (p ,) 

R 

where [(i + A + B + C) R] exp[j(tot - kR)) is the maximum possible 
pressure amplitude and the distances D 2 ^ D 3 ^ and 

measured from the origins to the particular element. The 
directivity is; 

U(G,<I>) = {[cos(kD,sin0sin<I)) + Acos(kD,sin0sin<l> - (p,) (14) 

(1+A+B + C) * * 

+ Bcos(kD^sin0cos<I> - (p^) + Ccos(kD_jSinGcos<I> + 94 )]“ 

+ lsin(kD|Sin0!>in<I>) - Asin(kD,sin0sin<I> - <p,) 

- Bsin(kD^sinOcos<I> - <p^) + Csin(kI)^sinGcOb<P + <P 4 ) 1 *} 
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z 




Figure 2.7 Quadrupole Geometry 

If D]^ = D 2 = D 3 = D 4 then comparison of Figs. (2.7) and 

(2.6) shows that the configuration is simply the Tesseral 
quadrupole of Fig. (2.6) where has been shifted 45*^. 

Additionally, the axial quadrupole of Fig. (2.5) is 
obtained using Park's program by placing source B at an 
infinitesimally small distance from the array center 
(i.e., d 3 = 0.00001 meters for any A ) with amplitude 

weighting equal to twice sources 1 and A and source C 
amplitude weighted to zero, which removes its contribution 
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from the array. This axial quadrupole can also be obtained 
by placing sources B and C at the array center and amplitude 
weighting equally. 

These equations were utilized to provide verification 
data for the laboratory experiments. The three element 
array data were verified by amplitude weighting element C of 
Eq. (14) to zero. Dipole data were verified by amplitude 
weighting elements B and C of Eq. (14) to zero. 
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III. RXPRRTMENTAT. PROCEDURE AND DEVELOPMENT 



A. BEAM PATTERN MEASURING SYSTEM OVERVIEW 

The beam pattern measuring system is an instrumentation 
package which was designed to power and rotate an acoustic 
array and then receive the signal, compute the pressure 
distribution of the array and display this data in both a 
tabulated format and plot of sound pressure level in 
decibels versus degrees of rotation. All data were taken in 
an anechoic chamber to minimize surface interference. 

The system used for measuring beam patterns of different 
arrays is divided into two major groups. These are the 
signal generation and input section and the automated 
advancement, data collection and display section. An 
overview of these is now presented along with a block 
diagram of the system. Fig. (3.1). A detailed description 
of each piece of equipment is contained in Appendix A. 

The Hewlett-Packard (HP) Function Generator Model 3314A 
was utilized for all signal generation. The signal passed 
through an individual HP-467A Power Amplifier to each of the 
University Sound Model ID-30 acoustic driver units. Setting 
the output levels of the driver units was accomplished 
utilizing the variable gain control on the power amplifiers. 
This constituted the signal generation and input section. 
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Figure 3.1 System Diagram 
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The receiving and processing section was controlled by 

the HP-86B computer. It directly controlled the array 

rotation and the reading of all measurements of peripheral 

instruments via the Hewlett-Packard Interface Bus (HPIB) . 

The computer sent clockwise stepping commands to the motor- 

power/direction controller via the HP-3421A Data 

Acquisition/Control (DAC) to rotate the array in stepped 

increments of about one degree. Following each advancement 

of the array, the computer read both the resistance output 

of the array rotation measurement potentiometer and the 

voltage output of the General Radio Model 1962 Microphone. 

The array rotation potentiometer was a 0 to 10 Rohm 

potentiometer mounted on top of the rotation motor and 

measured the angle of rotation of the array (0 to 10 Kohms 

corresponded to 0° to 360°) . These resistance and voltage 

readings were then stored in the computer for use in the 

data analysis and plotting routines. Once the data was 

o 

taken for 360 of rotation, it was used to compute the 
rotation angle and its corresponding decibel (dB) level for 
each datum. These values were then output to the computer 

monitor and plotted in the form of a polar plot of dB versus 

degrees of rotation. 

B. BEAM PATTERN PROGRAM 

The Beam Pattern program listed in Appendix B, consists 
of a series of sub routines. The program starts by 
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providing addresses for the two multimeters and the data 
acquisition/control unit. Next the initialization 
subroutine resets and initializes these three instruments. 

Once the equipment is ready to operate, the incrementing 
subroutine is called. This is an iterative and interactive 
routine used to rotate the array and measure the voltage of 
the sensed pressure field and the corresponding resistance 
reading which is proportional to the angular position of the 
array. Up to 360 data points can be collected per data run. 
After the initial datum is read, the array rotator motor is 
energized for 0.5 milliseconds to increment the array 
azimuth. The array is then allowed to coast down for 4 
seconds prior to reading the voltage and resistance. This 
selection of rotate and coast down times is discussed later 
in the chapter. Once the resistance and voltage readings 
are taken, they are stored in an array for later transfer to 
a data disk and for use in the plot subroutine. 

Finally the plot subroutine calculates the relative 
decibel level of the beam pattern based on the measured 
voltages. It then takes the resistance reading and converts 
it to the angle of rotation. The calculated decibel levels 
are then normalized to zero decibels based on the highest 
received voltage reading. Finally, the data points are 
plotted on a polar plot of relative sound pressure level in 
decibel versus angle of rotation. 
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All raw data is stored on a data disk after a complete 
data run. Each data run is stored by run number, as shown 
in Table II. It can be retrieved and printed at the 
printer as received voltage and resistance for each datum 
point or printed and plotted as decibel versus angle of 
rotation . 

The physical weight of an array was a major variable in 
this program with regard to the selection of the number of 
data points needed. Arrays weighing less than 4 pounds 
required 360 data points for a full 360^ of rotation. The 
number of data points required for larger arrays was 
decreased as array weight increased. For the four element 
array, weight was substantial enough (19 pounds) to require 
only 317 data points for a full 360*^ rotation. This was 
caused by the fact that the heavier arrays would coast-down 
for a larger distance after the 0.5 millisecond motor-on 
period . 

Another contributor to variation in the required number 
of data points was array balance. When utilizing a three 
element array, the array was sufficiently out of balance to 
require a greater number of data points for full rotation. 
This was attributed to minor binding of the rotator shaft as 
the array coasted down. 

Selection of motor run time and coast-down time was 
based on experimental investigation. A 0.5 millisecond 
motor on time provided 1 rotation regardless of coast-down 
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time with a light weight array. The 4 second coast-down was 
selected because smaller periods did not allow the array to 
stop and the resistance multimeter to stabilize prior to it 
being read. Motion at the variable potentiometer when the 
resistance values were read yielded erratic reading on the 
multimeter due to electromagnetic interference generated by 
the motor. 

C. EXPERIMENTAL APPARATUS DEVELOPMENT 

The experimental apparatus was designed to optimize a 
thorough search of two, three and four element arrays and to 
establish their respective beam patterns. Previously, 
several experimental investigations have been conducted 
utilizing a two microphone array with a single source. It 
was, therefore, decided that reciprocity would be invoked 
and the array would consist of multiple ID-30 driver units 
(Appendix A) with a single microphone utilized for sensing 
the pressure field. 

Six ID-30 driver units were checked individually in the 
anechoic chamber for their beam pattern impedence, 
resistance and reactance characteristics. This data was 
utilized to select the best matched set of four driver 
units. All drivers were examined between 5 Hz and 3000 Hz 
in an effort to locate and identify the resonances in the 
individual speakers. Fig. (3.2) shows a representative 
resistance curve for these drivers between the above 
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Figure 3.2 Typical ID-30 Resistance Plot 



frequencies. It was desired to conduct all analysis at 
frequencies in relatively flat areas away from resonances. 
This was done to assure phase stability since reactance is 
not varying. The speakers utilized in this investigation 
had resonances as per Table I. 

TABLE I 

ID-30 DRIVER UNIT RESONANT FREQUENCIES 
Speaker Number ; First Resonance (Hz) ; Second Resonance (Hz) 



1 


689 


2290 


2 


660 


2225 


3 


780 


2150 


4 


798 


2129 


5 


740 


2253 


6 


778 


2125 



Matched set number 1 consisted of speakers 1 and 2 and set 
number 2 of speakers 3 and 4. All array beam patterns were 
run at 300 Hz and 960 Hz which fell in the flat areas of the 
impedance curves. 

The array existed in two individual configurations, a 
linear and a crossed axis system with variable tilt (Fig. 
(3.3)). The linear system was used for investigation of 
doublets, some triplets and quadruplets (in the form of an 
axial quadrupole) . The crossed axis array proved to be 
extremely versatile and useful, providing four arms, each 
70 centimeters in length, for holding the drivers. 



Additionally, it provided a means to vary the angle© a full 
90°. Each driver was equipped with a harness assembly which 
mounted to the array rod in such a manner that the driver 
face was located 2.5 centimeters off the axis. This did not 
prove to be a matter of concern since all pressure 
measurements were made in the far field. 




Figure 3.3 Array Configurations 
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A scattering and shading problem of small magnitude was 
encountered because of the large center piece on the crossed 
axis array. This scattering problem was largely resolved by 
wrapping the center piece with fiberglass. 

The General Radio Model 1962 microphone was located 3.2 
meters from the center of the array placing it in the far 
field. It was mounted at the same height in the anechoic 
chamber as the array center. Changes in tilt angle theta 
were accomplished by tilting the array. 

The outputs of the ID-30 driver units were balanced 
prior to each data run. It was done by aligning the array 
in the vertical plane with the drivers facing the 
microphone. Then the output of each driver was adjusted 
utilizing the variable gain control on the power amplifier. 
This allowed all the drivers to be driven within a received 
voltage of 0.002 VAC of each other out of a typical level of 
1 Volt. 

Output signals from the HP-467A Power Amplifiers were 
monitored on an oscilloscope as was the received signal from 
the microphone. Because the signals were clean and without 
any substantial noise, filters were not utilized in the 
microphone circuit. 
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IV. DATA ANALYSIS AND VERIFICATION 



After establishing the ability to accurately rotate an 
acoustic array while collecting and displaying sensed 
pressure levels under computer control, it was necessary to 
ensure the data were accurate. This could be done by taking 
a set of parameters and obtaining both theoretical and 
experimental beam patterns for them. Verification was 
conducted utilizing the results of a FORTRAN program 
developed by Park (Ref. 1) which yields a variety of useful 
plots with the help of DISSPLA. 

A. THEORETICAL BEAM PATTERN PREDICTION 

With the program developed by Park, quick beam pattern 
plots could be made to compare with experimental patterns. 
As an example. Fig. (4.1) depicts a three dimensional 
representation of a dipole. Clearly, this particular is not 
of much use. 

The three dimensional plot was sectored for greater 
utility. Fig. (4.2) shows a sector representation for a 
four element array with kD^ = kD 2 = 3.52 and kD^ = kD^ = 0. 
Note that the kDj^ utilized in the remainder of the thesis 
related directly to the spacing from the center of the array 
to the element as related to Fig. (2.7). Refer to Fig. 
(2.7). The Z-axis of the plot is beam pattern (dB) , the Y- 
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DIRECTIVITY PATTERN 




kDi=kD2= 3.52 
kD3=kD4= 0 

Dl=02 D2=02 
D3=0. D4=0. 

VIEW ANGIJ:= 20, 20 
P2=180 P3=0 P4= 0 
A=l. B=0. C=0. 



Figure 4 . 1 Three Dimensional Dipole Representation 
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Figure 4. 




kD2=kD2=kD3=kD^= 10.55 

pi=0.6 D2=0.6 
03=0.6 D4=0.^ 

^EW ANGLE=+30.+60 
P2= 0 P3= 0 P4= 0 
A=l. B=l. C=l. 



2 Sectored Three Dimensional Plot ( 4 Element ) 
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axis psi (degrees) , and the X-axis theta (degrees) . 
Directivity is scaled from 20 to -40 dB (with 0 dB being the 
highest achievable value because of normalization) , psi from 
180*^ to 360°, and theta from 0° to -90°. This plot appears 
to be very busy but proves to be quite useful. It is easy 
to evaluate this sectioned representation by following the 
contour line of the directivity above the line theta equals 
0 ^ for psi from 180° to 360°. This provides the same dB 
levels at 180° < psi < 360° as the polar plot representation 
plotted in Fig. (4.3). 

B. EXPERIMENTAL DATA-SINGLE ELEMENT DIRECTIONALITY 

As Stated in Chapter III, testing was conducted at 300 
Hz and 960 Hz which were in the flat areas of the drivers' 
resistance curves, away from their resonance peaks. This 
was done to assure phase stability between drivers since the 
reactance is not varying (large phase shifts do not occur 
except near resonance) . 

Initial data runs were conducted with a single ID-30 
acoustic driver to establish the directionality of each of 
the drivers. The drivers were rotated 360° at three 
separate tilt angles of the driver face, © = 0°, 45° and 90° 
A representative graph of the directionality is contained 
in Fig. (4.4) which was made with the driver at a tilt angle 
of 0° and pointing at the microphone prior to rotation. The 
dotted line is the laboratory beam pattern. Note that there 
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PHI VS. DIRECTIVITY 



0 




THCTA= 0 DEGR^^ 

kD]^=kD2=kr)3=l^-4" 



Ul=0.6 D2=0.6 
D3=0.6 D4=0.6 
P2=0 P3=0 P4= 
A=l. B=l. C=l. 



10.55 

0 



Figure 4.3 Polar Plot of a Four Element Array 
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Figure 4.4 ID-30 Acoustic Driver Directionality 



is a 1 to 2 dB drop below the 0 dB level as the driver 
points away from the receiver which shows the driver is not 
completely omnidirectional. It should also be noted that 
the program's normalization routine sets the highest 
received microphone voltage level to 0 dB and normalizes all 
other voltage levels in the beam pattern to this level. 
This was done to provide greater ease in plotting of beam 
patterns so as to produce a uniformity between plots. 

C. EXPERIMENTAL AND THEORETICAL DATA - VARIOUS ARRAY 

CONFIGURATIONS 

An investigation of two, three and four element arrays 
was conducted. For each of these configurations, several 
different combinations of spacing and amplitude and phase 
shading were looked at. Table 2 at the end of this chapter, 
is a list of these variations for which beam patterns were 
obtained. Note that phase shading was limited to either § 
or 180®. In viewing the figures, the term "x-pole" is used 
loosely to represent the term "x-let" (i.e., quadrupole for 
quadruplet) . 

For the two element array, beam patterns were run in the 
laboratory for both- in- phase and 180® out-of-phase doublets. 
An example of the out-of-phase doublet beam pattern obtained 
in the anechoic chamber is shown in Fig. (4.5). (Note that 
in this figure kd is the quantity 2kDj when Dj^ =02-) 

Its theoretical verification, obtained from Park's program 
is shown in Fig. (4.6). A contributing factor to their 
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DIPOLE BEAM PATTERN (dB) 
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THCTA VS. DIRECTIVITY 

0 




PHI= 90 DEGREE 

kDi=kD2= 5.276 , kD3=kD4= 0 

D1=D2= 0.5 
D3=D4= 0 

P2=180 P3- 0 P4= 0 
A=l. D---0. C=0. 



Figure 4.6 Theoretical Dipole Plot 
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small differences, is the slight directionality associated 
with the acoustic driver. The slight offset of the 
laboratory plot nulls is caused by the array not being 
exactly aligned to broadside at 0 Kohm on the potentiometer. 
The 90*^ rotation between the two plots is caused when theta 
is fixed at 90® vice 0® for theoretical plots. Because of 
the narrowness of the nulls, laboratory data would require a 
sample exactly with the array broadside to achieve the null. 
Other than the difference in nulls, there is exceptionally 
good agreement between experimental and theoretical beam 
patterns . 

In the investigation of three element arrays, the arrays 
were configured as equilateral triangles, isosceles 
triangles and as a linear array. For phase shading, three 
cases were looked at: all elements in phase, the center 

element 180 out- of-phase, and one end element placed 180 
out-of- phase . These arrays were also subjected to various 
amplitude shading (the axial quadruplet composed of three 
elements with the center element amplitude shaded to twice 
the value of each of the end elements will be discussed in 
the quadruplet section) . Experimental data for a three 
element isosceles triangle with all elements in phase and 
= kD 2 = 2kD3 = 5.28 are shown in Fig. (4.7) with its 
theoretical verification shown in Fig. (4.8). The deep 
theoretical nulls seen on Fig. (4.8) are not present in the 
laboratory data. This can be attributed to a small variance 
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TRIPOLE BEAM PATTERN (dB) 
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THETA VS. DIRECTIVITY 




PHI= 90 DEGREE 
kDi=kD2-2kD3= 5.276 

DI-03 D2=0.3 
D3=0 15 D4=0. 

P2= 0 P3= 0 fH=0 
A=l. B=l. 0=0. 



Figure 4.8 Three Element Isosceles Plot 
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in the equal amplitude shading, a small amount of scatter 
from the large array center joint, and a slight 

directionality of the individual acoustic drivers. 

The beam patterns of a triplet in the form of an 
isosceles triangle are shown in Figs. (4.9), (4.10), and 

(4.11) to show the effects of tilting the array out of the 

horizontal to angles of theta equal to 30*^ , 60*^ and 9^ 

respectively. The plane of the elements was initially 

perpendicular to the x-y plane, making theta, the angle 

o 

between the z-axis and the array plane equal to 0 . The 

array plane was then tilted 30^ about the x-axis and rotated 
360° about the z-axis to obtain the data for Fig. (4.9). 
The array plane was then tilted in 30 increments to obtain 

Figs. (4.10) and (4.11). For this triplet case, kD ^ = kD ^ = 

2kD3 = 5.28 with kD3 sensing the elevation change (Fig. 

(4.12) ) . 
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TRIPOLE BEAM PATTERN (dB) 
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TRIPOLE BEAM PATTERN (dB) 
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Figure 4.11 Beam Pattern for 90 Tilt in Theta 





Figure 4.12 Array Configuration for Tilt Theta 

These variations in tilt angle (theta) show that the 3^ 
tilt has a diminished dB level for all lobes except the 
broad lobe coming into and exiting broadside. This is due, 
in large part, to the directionality of the drivers, because 
they are facing the microphone to a greater extent at 3^ 
than at 90° of tilt and mutual shading exists. 
Additionally, for these various tilt angles, the Z-axis is 
rotating about the local vertical, as explained in section D 
of this chapter. 

The last major configuration investigated was the 
quadruplet. It was looked at both as an axial quadruplet. 
Fig. (2.5), and as a planar array. Fig. (2.7). An 
experimental beam pattern for a planar array quadruplet with 

*^^1 = kD2 = kD3 = kD4 = 5.26 and all elements in phase is 
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shown in Fig. (4.13) along with a theoretical beam pattern, 
Fig. (4.14). There is good agreement between the two plots. 

A plot of an axial quadruplet represented by an 
amplitude shaded triplet, obtained in the laboratory is 
shown in Fig. (4.15). Comparing Fig. (4.15) to Fig. (4.13), 
it is seen that when two elements of the planar array are in 
line, the arrival from the more distant source is weakened 
by the shadowing effect of the nearer source. Therefore, the 
near source can be emphasised by as much as one dB when 
compared to the far source. Additional plots of various 
configurations are contained in Appendix C, where good 
agreement continues to be seen. 

D. EFFECTS OF ARRAY TILT, THETA 

Referring back to Fig. (4.2), it is again noted that 
contour lines are directly related to polar plots when 
either theta or psi is fixed. However, the countour 
representation is not as useful when theta and psi are both 
changing with rotation of the array as in the case where the 
array has been tilted in theta. In this case the polar 
curve is found by superimposing a trajectory which is the 
specific locus of© and ^ in the © plane onto the 
contoured surface. 

For this work, the microphone was fixed in the far field 
at a distance of 3.2 meters and level with the center joint 
of the array. Tilting the array in angle theta caused an 
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THETA VS. DIRECTIVITY 



0 




PHI= 90 DEGREE 

kDi=kD2=>^3=^D4= 5.215 

Dl=05 D2=05 
D3=0.5 D4=05 
P2=0 P3=0 P4=0 
A=l. B=0. C=0. 



Figure 4.14 Theoretical Planar Quadruplet Plot 



55 



TRIPOLE BEAM PATTERN (dB) 



I 







56 



R-3.20 M kDi=kD2=kD3= 4.713 
PHASE-0.100.0 Deg AMPL- 0.491 Vac SHADE- 1. 2. 1 



offset between the Z-axis and the local vertical about which 



the array rotated. This can be seen in Fig. (4.16). In the 
left sketch, the array elements are aligned in the standard 
rectangular coordinate system, with the elements in the x-y 
plane and rotating about the z-axis. However, when the 
plane of the array is tilted 30® about the x-axis, Y' = Y+3^* 
and z ' = z+3cP. A 90® rotation about the z-axis to 27^ 
yields z ' = 90® - 30® = 60® . This means that a 36^ 
rotation about the z-axis causes Yl and Zl to vary from 12^ 
to 60® and back to 120® . This is due to the fact that the 
microphone is fixed in space in the far-field rather than 
tilting with the plane of the array. 





Figure 4.16 Z-axis Offset Effects 

This problem does not occur in the classical representation 
where the microphone is fixed with respect to the axis of 
rotation of the array. 

E. ERRORS 

The sources of errors in the beam pattern measuring 
system were examined both mathematically and empirically. 
Sources of error for direct error analysis were based on 
multimeter values and power amplifier settings. With output 
reading from the multimeters and power amplifier exhibiting 
precision to the third decimal place, instrumentation errors 
from these sources are less than 1% and are not major 
contributing errors to measured beam patterns. 
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The major source of error in the system was determined 

empirically to be the acoustic driver positioning on the 

array. The effects of a small change in kd are shown in 

Figs. (4.17), (4.18), and (4.19). Here, it is seen that 

small changes in kD^ prove to have a large effect on the 

size of the small lobe at ± 90*^ and 0° and 180° , thus 

yielding the largest error. A change of ± 5% in kd effected 

the structure by up to 6 dB as seen from the behavior of the 

o 

small lobe at ±_ 90 . This change relates to a positioning 
error of about ± 3 centimeters and presented a small problem 
with respect to exact positioning of the speakers. 

F. VALIDATION OF RESULTS 

The laboratory beam patterns prove to be in good 
agreement with theoretical beam patterns. Small uniform 
beam pattern amplitude differences are merely a consequence 
of the difference in normalization methods employed in the 
theoretical and experimental beam pattern processing. The 
theoretical program normalized the directivity to the sum of 
the amplitude weights, whereas the laboratory program 
normalized the beam pattern to the highest received voltage. 

An independent verification of the experimental program 
was performed by a thesis group under the auspices of 
Professor Don Walters of the Naval Postgraduate School 
Physics Department, who asked to use this experimental 
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THETA VS. DIRECTIVITY 



0 




PH1= 60 DEGREE 

kD]^=kD2=kD3=kD4= 4.75 

Dl=0.5 D2=05 
D3=0,5 D4=0.5 
P2= 0 P3= 0 P4= 0 
A=l. B=l. C=l. 



Figure 4.17 kD^= 4.75 for Error Corputation 
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THETA VS. DIRECTIVITY 

0 




PHI= 60 DEX5REE 
kDj=kD2=^3=1^4= 
U1=0 5 D2=0.5 
D3=0.5 D4=0.5 
P2= 0 P3= 0 P4= 0 
A=l. B=l. C=l. 



Figure 4.18 kDj^- 5.0 for Error Conputation 
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THETA VS. DIRECnVTTY 



0 




PHI= 60 DEGREE 

kDi=kD2=kD3=kD4= 5.25 

Dl=05 D2=05 
D3=05 D4=0.5 
P2= 0 P3= 0 P4= 0 
A=l. B=l. C=l. 



Figure 4.19 kDi= 5.25 for Error Conputation 
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The laboratory system was placed 



system 
in use 
for a 5 
further 



for related work. 

by this research group in the selection of drivers 
X 5 array currently under investigation and will be 
utilized to establish beam patterns for this array. 



TABLE II 



ARRAY CONFIGURATIONS 



Run No/ 










Amplitude Spacing 


No. 


Freq 




Phase Tilt 


Shad ing 


from Center 


Elements ; 


: (Bz) 


I 


(deg) I (deg) i 


(VAC) 


1 (CM) 


1 


300 






0 






1 


300 






45 






1 


300 






90 






1 


960 






0 






1 


960 






45 






1 


960 






90 






2 


300 




0,0 


0 




30 


2 


300 




0,180 


0 




30 


2 


960 




0,0 


0 




30 


2 


960 




0,180 


0 




30 


2 


960 




0,0 


90 




30 


2 


960 




0,180 


90 




30 


2 


300 




0,0 


90 




30 


2 


300 




0,180 


90 




30 


2 


300 




0,0 


45 




30 


2 


300 




0,180 


45 




30 


2 


960 




0,0 


45 




30 


2 


960 




0,180 


45 




30 


4 


960 




0,0, 0,0 


0 




30 


4 


300 




0,0, 0,0 


0 




30 


4 


960 


0 


,180,3,180 


0 




30 


4 


300 


0 


,180,0,180 


0 




30 


4 


960 




0,0, 0,0 


30 




30 


4 


300 




0,0, 0,0 


30 




30 


Quadll/4 


300 


0 


,180,0,180 


30 




30 


Quadl2/4 


300 


0 


,180,0,180 


60 




30 


Quadl3/4 


960 


0 


,180,0,180 


60 




30 


Quadl4/4 


960 




0,0, 0,0 


60 




30 


Quadl5/4 


300 




0, 0,0,0 


60 




30 


Quadl6/4 


960 




0,0, 0,0 


90 




30 


Quadl7/4 


960 


0, 


,180,0,180 


90 




30 


Quadl8/4 


960 


0, 


,180,0,0 


0 




30,30,30,0 


Tr i 1/3 


960 




0,0,0 


0 




30 


Tri2/3 


300 




0,0,0 


0 




30 


Tri3/3 


300 




0,180,0 


0 




30 


Tri4/3 


960 




0,180,0 


0 




30 


Tri5/3 


960 




0,0,180 


0 




30 


Tr i6/3 


960 




0,0,180 


30 




30 


Tr:-^/3 


960 




0,180,0 


30 




30 


Tri8/3 


960 




0,0,0 


30 




30 


Tri9/3 


960 




0,0,0 


30 




30,15,30 


Tril0/3 


960 




0,0,0 


0 




30.15.30 


Trill/3 


960 




0,180,0 


0 




30,15,30 


Tril2/3 


960 




0,180,0 


30 




30,15,30 


Tril31/3 


960 




0,180,0 


60 




30,15,30 


Tril4/3 


960 




0,0,0 


60 




30,15,30 


Tril5/3 


300 




0,0,0 


60 




30,15,30 


Tril6/3 


960 




0,0,0 


90 




30,15,30 


Tril7/3 


960 




0,180,0 


90 




30,15,30 


Tril8/3 


960 




0,0,180 


90 




30,15,30 


Tri20/3 


960 




0,0,0 


0 


1,1,1 


30,0,30 


Tri21/3 


960 




0,180,0 


0 


1,1,1 


30,3,30 


Tri22/3 


960 




0,0,0 


0 


1,2,1 


30,0,30 


Tri23/3 


960 




0,180,0 


0 


1,2,1 


30,0,30 


Tri24/3 


960 




0,0,180 


0 


1,2,1 


30,0,30 


Tri25/3 


960 




0,0,0 


0 


1,2,1 


30,30,30 


Tri26/3 


960 




0,0,0 


30 


1,2,1 


30,30,30 


Tri27/3 


960 




0,0,0 


60 


1,2,1 


30,30,30 


Tri28/3 


960 




0,0,0 


90 


1,2,1 


30,30,30 


Tri30/3 


960 




0,0,0 


90 


1,4,1 


30,30,30 


Tri31/3 


960 




0,180,0 


90 


1,4,1 


30,30,30 


Tri32/3 


960 




0,0,0 


0 


1,2,1 


20,20,20 


Tri33/3 


960 




0,180,0 


0 


1,2,1 


20,20,20 


Tri34/3 


960 




0,180,0 


0 


1,2,1 


20,10,20 


Tri36/3 


960 




0,0,0 


0 


1,2,1 


20,10,20 


Tri37/3 


960 




0,0,0 


0 


1,2,1 


10,10,13 


Tri38/3 


960 




0,180,0 


0 


1,2,1 


10,10,10 


Tri39/3 


960 




0,180,0 


90 


1,1,1 


20,30,20 


Tri40/3 


963 




0,180,0 


90 


1,2,1 


17.9,17.9,17.9 


Tri41/3 


960 




0,0,0 


90 


1,2,1 


17.9,17.9,17.9 


Tri43/3 


960 




0,0,0 


90 


1,2,1 


26.8,26.8,26.8 


Tri44/3 


960 




0,180,0 


90 


1,2,1 


26.8,26.8,26.8 



Note: (1) Additional data runs were made in some 

of the same configurations with fiberglass 
wrapped around the centerpiece to identify 
the effects of scattering. 



(2) Additional data runs were made in some 
of the same configurations to verify 
repeatability. 



64 



V. CONCLUSIONS RECOMMENDATIONS 



The doublet forms the basis for countless sonar and 
radar arrays currently utilized in both military and 
civilian applications. It is the ability to control the 
amplitudes and phases of the element of an array which 
allows the tailoring of the array radiation pattern to 
achieve the desired element of interest and promotes the 
maximization of array performance by increasing signal, 
decreasing noise or both. 

The ability to rapidly measure the beam pattern of an 
array with complete confidence in the results is the basis 
behind this study. A laboratory system was developed to 
allow collection of a 360 point beam pattern in a 45 to 50 
minute period. Software is easily modified to allow for 
collection of more or less than 360 data points. One of the 
major advantages of the computer-controlled method of 
obtaining a beam pattern is the relatively small amount of 
time required to achieve an accurate beam pattern based on a 
large data base. 

Recommendations for improving the current laboratory 
system are based on reducing electrical noise of the 
rotation motor and reducing the size and weight of the 
array. First, the noise of the rotation motor prohibits 
continuous rotation while reading the resistance on the 



65 



potentiometer. Since a CW signal does not require broadband 
detection, filtering could be accomplished by use of a lock- 
in, fast fourier transform or heterodyne filter. Second, the 
array weight should be reduced to lessen the drag on the 
rotation motor and to reduce the amount of coast-down time 
required for heavier arrays. Next, the array center joint 
should be reduced in size to lessen any shadowing effects 
from the large, metal center joint. (Scattering effects 
seemed to be small judging from the lack of any effect on 
the polar plots when this joint was wrapped in fiberglass.) 
Also, driver mounts need to be designed so that the center 
of the driver face is directly on the axis. The 2 
centimeter offset utilized for this investigation was 
minimized by making all readings in the far field. Last, 
greater distance between the driver and microphone would 
also help to lessen the 1 dB variation seen when an array 
was in an endfire orientation. This only affected the depth 
of the nulls. 



66 



APPENDIX A 



EQUIPMENT DESCRIPTION 

1. ANECHOIC CHAMBER 

All experimentation was conducted in the anechoic 
chamber, Room 019, Spanagel Hall, Naval Postgraduate School. 
The anechoic chamber is a 99% echo free enclosure for sounds 
above 100 Hertz. Acoustics research can be conducted on 
sound sources, receivers and scatterers with minimal 
interference from surface reflections or external noises. 
The size of the outer room in which the anechoic chamber 
floats is 37 feet x 24 feet x 21 feet high. Noise 
transmission from outside is minimized by the floating room 
within a room construction. The outer 12 inch concrete- 
walled room is separated from the inner room of concrete- 

block sides and floor by a, 2 inch blanket of fiberglass or 

cork. Forty inch fiberglass wedges are attached to all 
surfaces of the inner room. These are designed to trap and 

absorb sounds which are incident upon them. A suspended 

floor is constructed of a grid of 225 wire cables, each 
stretched at a tension of 150 to 200 pounds. 

2. HP-86B PERSONAL COMPUTER 

The HP-86B computer, operating in Hewlett-Packard BASIC 
language, is equipped with a 64 Kilobyte (K-byte) memory 
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with two additional 128 K-byte memory modules. The system 
proved more than sufficient for both program development and 
as an instrument controller. Instrument control by means of 
the computer was accomplished with the HP-IB (Hewlett- 
Packard Interface Bus) which is a standard connection on 
most HP equipment. This facilitated instrument control for 
the various system components. 

3. HP 7470A GRAPHICS PLOTTER 

The HP 7470 a Graphics Plotter is a vector plotter which 
utilizes HP's micro-grip drive technology, moving both the 
paper and the pen. Multicolor plotting is available, but 
there are only two pen positions so pauses must be used to 
change pen colors. Only one line of program code is 
required to shift the output to and from the plotter to the 
HP-86B CRT. 

4. HP-4192A LOW FREQUENCY IMPEDANCE ANALYZER 

The HP-4192A LF Impedance Analyzer is a fully automated 
instrument designed to measure a large range of impedance 
related characteristics. It can be used to measure eleven 
impedance parameters; absolute value of impedance Z, 
absolute value of admittance Y, phase angle©, resistance R, 
reactance X, conductance G, susceptance B, inductance L, 
capacitance C, dissipation factor D, and quality factor Q. 
The HP-4192A was used to measure the frequence response. 



68 



impedance, resistance and reactance of the University Sound 
ID-30 driver units as discussed in section 10 of this 
appendix . 

5. HP-3421A DATA ACQUISITION/CONTROL UNIT 

The HP-3421A Data Acquisition/Control Unit (DAC) has a 
5.5 digit voltmeter, a 10 kHz counter and an HP-IB 
controller interface. The DAC was utilized as a computer 
controlled switching unit with channels 1 and 2 configured 
as actuators. This provided a voltage path to control 
relays in the motor power and direction controller for 
rotation of the ID-30 speaker array. Fig. (A.l). 

6. HP-3478A MULTIMETER 

The HP-3478A Multimeter provides both resistance and 
voltage measurement capabilities with 3.5 to 5.5 digit 
resolution with + 0.3 mV accuracy on the 100 Hz to 20 kHz 
scale. The multimeter was remotely controlled by the 
computer via the HP-IB. Voltage readings from the GR-1962 
microphone and resistance readings from the variable 
potentiometer on the array rotator being read by the 
multimeter and then sent to the HP-86B computer via the HP- 
IB. 



7. HP-3314A FUNCTION GENERATOR 

The HP-3314A Function Generator is a multimode function 
generator which provides sine, square, and triangular 
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Data Acquisition / Control Unit 




Figure A.l. Data acquisition connection to motor-power 
controller. 
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functions as well as any desired arbitrary wave form from 
frequencies of 0.001 Hz to 19.999 MHz. It was operated in 
the free run mode providing a continuous 1 Vac sine wave to 
Hp-467A Amplifier's inputs. 

8. HP-465A AMPLIFIER 

The HP-465A Amplifier is a general purpose preamplifier 
and impedance converter from 10 megohms to 50 ohms. It has 
a selectable gain of 20 or 40 decibels and is stable over a 
frequency range of 5 Hz to 1 MHz. The amplifier was used to 
amplify the received signal from the GR-1962 microphone. 

9. HP-467A POWER AMPLIFIER 

The HP-467A Power Amplifier is a 10 watt peak power 
amplifier and a -20 to +20 void dc power supply. It has a 

wide band width with low dc drift from 5 Hz to 1 MHz. A 
variable gain control allows adjustments of gain between 
zero and 10 with resolution better than 0.1% of full power. 
Four HP-467A's, one for each ID-30 driver unit, were 
utilized in adjusting power inputs from the function 
generator . 

10. UNIVERSITY SOUND ID-30 DRIVER UNITS 

The University Sound ID-30 's are medium power horn 
driver units designed specifically for wide midrange 
response. Their tropicalized voice coils are 2 inches in 
diameter with linen base molded phenolic diaphragms with a 
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minimal 16 ohm impedance. The frequency response is 85 Hz 
to 7.5 Khz. They were used as the sound sources. 

11. GENERAL RADIO TYPE 1962 MICROPHONE 

The GR-1962 microphone is a 0.5 inch diameter electret 
condenser, microphone, providing a flat response at zero 
decibels between 20 Hz and 15 Khz. It experiences rapid 
drop off above 15 Khz and is nonrespons ive above 21 Khz. A 
General Radio Model 1560-P62 inline power supply was 
connected to the microphone 1560-P43. preamplifier to 
amplify the received signal. 

12. HP-6215A POWER SUPPLY 

The HP-6215A power supply provides zero to 30 volts dc 
or 0 to 500 Mamps . It was used to provide switching voltage 
for the Motor power relay in the direction controller. 

13. MOTOR POWER AND DIRECTION CONTROLLER 

The controller for the array rotator controlled power to 
the rotator and also the direction of rotation. With leads 
to the rotator assembly potentiometer, it also provided 
resistances which corresponded to a given angle (0 to 10 
Kohms translated into 0 to 360 degrees) . A +5 volt dc 
signal to IC-1 via channel 1 of the DAC provided on/off 
power to the rotator and +5 Volts to IC-2 via channel zero 
of the DAC provided clockwise rotation of the driver 
assembly unit (Figure A. 2). [Ref. 8] 
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ccw 




a) Circuit Diagram for Modified Motor Power and Direction Controller 




b) 360® / 1 0 kO Potentiometer c) 1 6 Pin DIP Relay 



Figure A. 2. Circuit diagram for modified motor-power/ direction 
controller [Ref. 9., p. 57], 
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APPENDIX B 



BEAM PATTERN COMPUTATION AND PLOT PROGRAM 



10 h:^[NT '■ BI.mM PhrTLRN COriPUl ATION PNU PLOP PROGPAh 

M FPim ” * 

1 : f 
13 I 

30 ppinr *t SET-UP 

J I I 



30 
4 0 
EO 
FjO 
6 1 
r: 

70 
00 
00 
100 
1 10 
1 1 1 
1 20 



PPiNF CHRi < 10 ) 

P'RlUr '’THIS PROGRAM CONTROLS THE FOLLGUING APPARATUS TO MEASURE, DISPLAr,' 
PFUNl "PLOl mND SluRE THE DIRECTIONAL CHARACTERISTICS OF .A SINGLET,” 

PR INI "DOUbl.ET, TRIPLET, AND OUnORUPLET ARRAY.” 



PRINT ' 
PRIM! ‘ 
PRINT ' 
PR IN I 
PRINT 

I 

DiSP ' 



INSTRUMENT ADDRESS 

HP. 3473 MULTIMETER (UOLTS) 724 

” HP^3470 MULT IMP lER ‘RES) 725 

” DnTh AQUISITIOLNCONTHOl.. 703 

INITIALIZE INSTRUMENTS" 



12 1 I 



150 REM preset hND INITIALIZE INSTRUMENTS*** 

140 OUTPUl /03 ;"RS" i RESET 

150 I'liiPUT 725 ;''P5R5T1" i FREE RUN 50K0HM5 

iEO OUTfMJT :''24 ;"F2RATr’ » FREE RUN AUTOPnNGE AC UOLTS 

151 ' 

170 LISP "INITIALIZE PROGRAM" 

171 I 

130 OPTION EASE 1 

1 U0 LLiM nlHN-.E^ 1 E0 ) , C ALC^ANGLE ( 1 50 ) , NAMES I 30 1 
200 

2 10 FF-INIFP IS 701 
220 INlFCrR I,'3,L,M 

2 50 Dirl 560 > ,UOLT ( 550 ) ,OHMS_DATA( 360 ) , X ( 560 ) , Yv 560 ) ,K( 560 ) 
2 40 M‘'iiLI-0 

2S0 MAZ.S STORAGE IS ".DmTAI" 

250 niZP "NUMBER OF DATA POINTS DESIRED'?" 

2/0 LLPP 10,500 
25 0 INf'UT M 

230 DISP "UnNl NEU DATA (Y/M)^" 

500 R5“EP 10,300 
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3II.T INPUT 04 

J20 If 04 = "Y'' THOM GOTO 550 
550 IF uf='M" IMFN GOTO 700 
540 GOTO 500 

550 PEN ' • » OEFmULT PnRP.TItTEFS **' 

351 I 

550 I MUaT INSERT FOR LINES 590, 400, 410 mNU 810 PRIOR TO START OF RUN 
5f;5 I 

590 CREATE "• INSERT NAME OF RUN *",3E0,1B 
400 hbSlbNfl I TO "* INSERT NmNE OF RUN •" 

410 NANEi-"* INSERT NAME OF FIlE »" i FILE NAME 
400 FLILI SrALE = 500 ! L0CK_II1 MUAC 
4 50 COf 10 = 580/ 1 0000 I CONUERl OHMS TO DEGREES 
440 5ET..UP: ' SET_UP MEASUREMENT PARnMETERS ••• 

450 PR INI CMIH <I0) 

4',!0 OISF "MANUALLY ENTER DESIRED FREQUENCY FOR SPEAKER DRIUER" 

500 DISP "ENTER 0 IF NO CHANGE" 

510 DISP "ENIER I TO CHANGE TIRE" 

520 BEEP I0,EO0 
530 INPUr C 1 CHOICE 
540 IF C=0 THEN GOTO 560 
550 IF i>l I HEN GOTO 750 
5E0 FRINT ChPi (10) 

570 OUIPUT 709 i"CLS0” ' ROIATE CW 

575 PR INI " I " , "0Hn5_DATA i. OHNb ) ” , "OOL.T ( VAC ) " 

5a0 FOR 1=1 TO M 

590 OUTPUT 709 R'CLSI" I POWER TO MOTOR ON 
GOO WmII .5 I MILLISECONDS 

B10 OUTPUT 709 ;"OPNI'' I POWER TO MOTOR OFF 
G20 WhIT 4000 

G30 ENTER 723 ; OHMS_DATA( I ) 

E40 ENIFR 724 ; UOLTd) 

050 PRINf I ,uHM5 .DhTA( I ) ,UOLT( I ) 

EiiO rjEAl 1 

bbO FOR f=l TO N 

700 PPRini 1 i OHMS_DATA( I ) ,UOL T(I ) 

7 I 0 NE r I 

720 hSSIGNtt 1 TO • 

730 OnsuB PLUT_ROUTINF 
74 0 GOTO 7 '70 

750 ENinf NhME'B i NEW TlTLEi 
7B0 Gu 1 u 4t'0 
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/70 LNU 

/uJ ri[£P HiHnr 15 DATA RUM DESIRED""' 

7GD DEEP ID. 500 

610 r55IGM« 1 TO "• INSERT MAME OF RUM 
S.:0 FOR 1-1 TO M 

piO PLr^Od 1 ; UHMS_ljArA( I ) ,UOLT( I ) 

610 ME/1 I 
650 r.'jSIGfU 1 10 

3E0 GOTO 'i’30 
db I * 
dCG ‘ 

670 PLO r_ROUT I ME : • PLOTS POLnR DATA 

87 1 t 
8 7: I 

860 PLOTTER 15 705 
630 GRAPH 
900 Cl L EAR 
910 f^Efi I 

920 S' nl E R0 ,G0 ,-B0 ,60 

930 yt^A IS 0 ,1 0 , -40 ^40 

94 0 V fi “ I S 0 , 1 0 , - 4 0 , 4 0 

950 CSIZE 3 

9bO LUiE ] vPE 1 

970 FOP L-0 TO 40 STEP 10 

980 Mu'^E 0 ,L 

990 t 'jr<G b 

1000 IrPEl. L-40 

lOiO fiEXT L 

10:0 l iriE TYPE 1 

1030 C8I3E 6 

10. 10 Mrji'E 0,50 

104 1 1 

IOT: I iriiERT LhBEl.S 1050. 1080 AMD 1100 PRIOR TO START OF RUN 
104 3 I 

1050 LABEL ;“OUADRUPGLE BEAM PATTERN (dB)" 

1060 CSIZE 3 



1 


li'i"L 0. 


45 








1 


1 ritO'l ; 


"H-3.:0 M" " ."f 


d -3 . 30" , 


" " ,"FREO = 300 H:" 






Mn"8 0. 


-50 








1 t 00 


irEFI ; 


f-HASE = 0 , 1 80 ,0 , 1 60 


Deg " , " 


" ,"nNFL= 0.093 Uac" 


" , " ANGLE = 30 DeQ 



1110 Ml'^'E 4 5.0 
1 1 :o LM-G 5 
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1 13l^ label ; "MIKE" 

I N'J LIME TYPE 4 

IIE0 FOR L=I0 TO 40 STEP 10 

1100 I'lOUE I ,0 

1170 FOR 1=0 TO 5E0 STEP 5 
1 160 DRAW L'O.uS ( I ) ,L*SIN ( I ) 

1 1 90 NEXT I 
1200 NEXT L 
1210 LIME TYPE I 
1220 PEN UP 
1230 FEN 2 

1240 FOR 1=1 TO N ' BE6IN5 BEAU PATTERN COMPUTATION 
1250 IF UOLTa) MOOLT THEN MUOL T = UOLT( I ) 

1260 NEXT I 

1270 FOR 1=1 TO N 

1290 K = 0HM5_DATA\ I )/ 10000*360 

1290 K1=K 

1300 ir- 011MS_.nATA( I )' 1 0 THEM K=KH-360/2S5 

1310 IF 0HN5_DATA(I )' 10500 THEN GOTO 1400 

1320 '.'ECIOR = 20'LbT ( UOl T( 1 )/MNOl.T ) U0 

1330 LIME TYPE 3 

13)0 IF 'JECTOR 1 THEN UECT0R = 0 

1350 PRINT I ,K , HECTOR 

1 360 [ )=UE0T0R'C09 (K ) 

1370 V( I ) = OECT0R«SIN (K) 

I 380 PLOT ' ( I > ,Y( r ) 

1390 PEN UP 
H0O NEXT I 
1410 RETURN 
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TABLE III 



SAMPLE PROGRAM DATA COLLECTION OUTPUT 



3 

4 

c 



OHhSJnTA (OHfiS) 
4 



36 

6 b 

79 



^.^OLT 
. .^ 4:14 
. 7G0Z3 
.753G: 
.70383 
.66367 



6 

7 

8 

9 

10 
1 1 
12 

13 

14 



106 


.61826 


123 


.57763 


1 SO 


.54058 


1 94 


.52056 


215 


.51957 


251 


.53613 


281 


. 567 j4 


308 


.60736 


330 


.65774 


369 


.72735 



16 

17 

18 

1 9 
20 

2 1 
-5 

n ~7 

24 

25 

26 

27 

28 

29 

30 



34 



401 


. 80624 


429 


.66951 


458 


. 93059 


438 


1 .00539 


522 


1 .07 1 L9 


543 


1 . 1 2759 


572 


1 . 1 69':7 


600 


1.244 


629 


1 .28825 


6-18 


1 . 3225 


cn 


1 .361'72 


704 


1 .40042 


726 


1 . 42822 


755 


1 . ]5G4 


786 


1 . 47833 


814 


1 .49626 


846 


1 . 50755 


86 9 


1.51197 


895 


1 .51 105 


92 4 


1 .50482 
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(..J L‘j c.-J rj 



TABLE IV 



SAMPLE CAIXXTLATED BEAM PATTERN OUTPUT 



I 


Degrees 


Decibels 


T 


1 . Jtj433398505 


3 1 . j 1 7 J 5 3 7 8 2 7 


J 


. 82 3 


34 . 25425^3403 




1 .296 


33 . 954 4 396329 


4 


2.015 


33.34b 1 423354 


5 


2.844 


32.6507973577 


5 


3.6)5 


32 . 232559509G 


7 


4.503 


31 .6428331572 


8 


5.7G 


31 .0727607067 


9 


B . 954 


30.738552392 


10 


7.74 


30.7220173563 


1 1 


9.03G 


30 . 9945357053 


1 2 


10.116 


31.4 86004604 3 


1 3 


1 1 .083 


32.0652076323 


] 4 


11.98 


32.77022 1 597 1 


15 


13.284 


33.6440053209 


1 G 


14.436 


34.5554233377 


1 7 


15.444 


35. 1946281369 


1 3 


16.498 


35.7871 038 1 67 


1 J 


17.568 


36.4558277425 


20 


1 S.792 


37 . 0088989b 1 G 


2 1 


19.54 9 


u ( . 4 b 2 1 b «L^ J 1 J 


J. 


20.592 


37.9183967977 



2 3 

J 

0 



34 

35 



2 I . b 
:: .503 
23.329 
24 .228 



39.3055441245 
33 .509 1 395402 
38.837050l-’lb 
39.09039284 



25.344 


39.3343026072 


26.136 


39.5050357255 


27.13 


33.674 7 4 9 9 ^ 7 


28.296 


39. 804564321 G 


29.304 


39.9092775356 


30.456 


39.97474536^2 


31 .284 


40 


*r TO 


35.9943ri56671 


. j j . 26 4 


39.9583276 1 
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RIPOLE BEAM PATTERN (dB) 



APPENDIX C 



Experijnental and Theoretical Plots 



n 

j 
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THETA VS. DIRECnVITY 




PHJ=90 






kDi=kD2=2kD3= 5.276 

D 1=0.3 03=0.3 
03=0.15 0-4=0. 
P2=1Q0P3=0 P4=0 
I A=l. B=l. C=0. 



Figure C.2 Theoretical Three Element Isosceles Plot 
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QUADRUPOLE BEAM PATTERN (cJB) 
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THETA VS. DIRECTIVITY 




PHI=0 DEG31EE 

kDi=kD2=kD3=--kD4= 5.276 

D1=05 D8=0 5 
D3=05 D4=05 

0 P3=180 P4=180 
A=l. B=l. C=l. 



Figure C.4 Theoretical Planar Quadrupole Plot 
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TRIPOLE BEAM PATTERN (dB) 
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TRIPOLE BEAM PATTERN (dB) 
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TRIPOLE BEAM PATTERN (dB) 
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QUADRUPOLE BEAM PATTERN (dB) 
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GUADRUPOLE BEAM PATTERN (dB) 
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